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Abstract 

An infinite particle system of independent jumping particles is con- 
sidered. Their constructions is recalled, further properties are derived, 
the relation with hierarchical equations, Poissonian analysis, and sec- 
ond quantization are discussed. The hydrodynamic limit for a general 
initial distribution satisfying a mixing condition is derived. The large 
time asymptotic is computed under an extra assumption. 
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1 Introduction 

Particle systems in the continuum are describing infinitely many particles 
located in the Euclidean space M. d . One may equip these systems with dif- 
ferent types of dynamics, deterministic as well as stochastic. First of all, we 
should mention the Hamiltonian dynamics and related problems concerning 
the derivation of kinetic equations for classical gases, see, e.g. [Sin89j . An- 
other dynamics strongly motivated by physical applications, is the gradient 
diffusion of infinitely many particles [Fri87j . |Spo86|. Note that in spite of 
serious efforts in both cases, the answers obtained till now are far from being 
complete. 

In order to identify further reasonable types of random evolutions in the 
continuum, we may look at the well established types of dynamics on lattice 
gas systems. There are two important classes of Markov dynamics, namely, 
Glauber and Kawasaki. Both are constructed in such a way that a given 
Gibbs measure (equilibrium state) on the lattice becomes an invariant mea- 
sure for the dynamics. The Glauber dynamics is a birth-and-death evolution 
of the lattice gas, contrary to the jump type evolution in the case of the 
Kawasaki dynamics. The latter stochastic dynamics are especially interest- 
ing for the study of hydrodynamic limits of interacting particle systems due 
to their a priori conservation law, the number of particles, and the resulting 
existence of a continuous family of invariant measures |DMP91j . Continuous 
versions of Glauber dynamics were introduced in |G1581j . [BCC02] . [KL05] . 
[KLR07J. Continuous analogs of the Kawasaki dynamics are random evo- 
lutions of particle systems in which individual particles jump in the space 
with rates leading to a Gibbs state in the continuum as an invariant measure 
|G1581j . [KLR07] . 

In this paper, we concentrate on stochastic dynamics where each particle 
performs a jump process independently of the others. This type of dynamics 
one might call independent jump process or free Kawasaki dynamics in the 
continuum. One of the aims of this paper is to demonstrate that already 
for the free Kawasaki dynamics the situation is technically non trivial and 
different from the gradient diffusions. Another aim is to lay a solid ground 
for future investigation of interacting systems. 

In order to better understand the chosen framework, one has to recall 
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the underlying motivation. Infinite particle systems are introduced as an ap- 
proximation for finite but very large systems. The underlying finite systems 
are systems of finitely many particles in a bounded subset (as, for example, 
a ball or a cube) of M. d . Each particle independently performs a Markov 
jump process. This underlying jump process we call in the following the one 
particle dynamics. Motivated by physics, one assumes that the particles are 
indistinguishable. Therefore, we describe the collection of positions of n par- 
ticles not by a tuple (x\, . . . , x n ) with x, G M d , but by the set {x\, . . . , x n } 
of points (excluding a priori coinciding positions) or by the integer valued 
measure S Xl + . . . + 8 Xn . These notations include the indistinguishability au- 
tomatically. We call this symmetrized collection of positions a configuration 
of particles. The latter interpretation, via measures, is called the empirical 
field. A key idea in the study of large particle systems is to substitute them 
by infinite volume systems in order to avoid boundary and finite size effects. 
The price to pay for this substitution is technical subtleties and difficulties, 
which we will describe in the sequel. 

Let us consider an initial configuration of positions for the infinite particle 
system, denoted by 7. In any finite observation window, 7 should look like 
a configuration of a large, but finite, particle system looked at in the same 
window. In the notation for configurations as sets this just means that for 
any open bounded subset A C M. d there can be only finitely many points in 
7 fl A. In the interpretation of configurations as non- negative integer valued 
measures, this means that 7 is a Radon measure. The construction of such a 
type of processes for general underlying one particle Markov process is given 
in |KLR09] and we shortly recall this construction in Subsection 13.11 In our 
case, each particle is equipped with an independent identically distributed 
exponential clock. If the clock of a particle rings the particle performs a 
random jump independent of the other particles. As one considers infinitely 
many particles, the construction is non trivial. First of all, the process cannot 
be started in any initial configuration. Actually, one may easily construct ini- 
tial configurations for which the process explode, in the sense that infinitely 
many particles appear in a finite volume. Secondly, the infinite volume pro- 
cess is not any longer a jump process in the classical sense, because in each 
time interval a.s. infinitely many jump events occur. An essential step in the 
construction is to show that in any finite time interval only finitely many 
jump events are visible in a bounded observation window, i.e., the number of 
particles in the window stays finite and only finitely many particles pass the 
window. In Subsection 13.11 we add to the results of |KLR09] the description 
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of the path-space measure corresponding to the process. One easily sees that 
all Poisson measures (Poisson random fields) with constant intensities are 
invariant measures with respect to the free Kawasaki dynamics. If the one 
particle jump rate is symmetric than these measures are also reversible. 

Till now we spoke only about processes starting in a single configuration 
of particles. Choosing an initial configuration at random w.r.t. a probability 
measure on the configuration space, one can derive further classes of pro- 
cesses. If, for example, one distributes the initial configurations w.r.t. an 
invariant measure, then one obtains a so-called equilibrium process. For the 
equilibrium processes w.r.t. reversible measures powerful methods of Dirich- 
let forms may be applied. In particular, existence can be shown even for 
general interacting systems, cf. [KLR07] . In the case of the free Kawasaki 
process, one can apply in addition second quantization techniques, cf. Sec- 
tion HJ] which give a full description of the L 2 -theory, also for non-symmetric 
jump rates. 

One speaks of a non-equilibrium process if the initial distribution is not an 
invariant measure. For the interacting non-equilibrium Kawasaki dynamics 
even the existence of such a Markov process is unknown. As was mentioned 
above, the free non-equilibrium Kawasaki dynamics exists for a large class of 
initial configurations. The aim of this paper is to go beyond mere existence. 
We want to study the ergodic properties of the process, i.e., we consider the 
large time asymptotic of the process in the sense of limiting distributions. 
As a next step, we consider the so-called hydrodynamic limit for the free 
Kawasaki dynamics. 

We shall expect that the difficulties to treat non-equilibrium processes de- 
pend on the class of initial distributions. Poisson measures for non-constant 
intensities form a class which seems to be the easiest to be handled in the case 
of the free Kawasaki dynamics. These distributions are no longer invariant 
measures, they are, in general, not equivalent to any invariant measure and 
there exists no semigroup theory for the free Kawasaki dynamics in L p -spaces 
with respect to these measures, cf. Remark 15.11 

Note that the one-dimensional distributions of the dynamics can be com- 
pletely described in terms of the one- dimensional distributions of the under- 
lying one particle jump process. Here a delicate moment is that the initial 
distribution of this one particle process is the intensity function for the start- 
ing Poisson measure. A natural class of intensities to consider should thus 
include intensities corresponding to the invariant measures, i.e., all constant 
intensities. In the infinite volume, this prevents us to work just with inte- 
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grable intensities which from a probabilistic point of view would be natural 
for the one particle process. A right framework seems to be the cone of all 
non-negative measurable bounded intensities, but it contains intensities for 
which there seems to be no asymptotic, cf. Remark 15 .81 The study of the 
large time asymptotic of the infinite particle processes reduces to the study 
of the large time asymptotic of the one-particle process for these initial in- 
tensities. 

The non-ergodicity of the infinite particle processes is reflected in the non 
ergodicity of the one particle processes in this class of initial intensities, see 
Subsection 15. II Under the additional assumption that the Fourier transform 
of the initial intensity is a signed measure we can show, applying Fourier 
analysis, that all processes with Poissonian initial distribution with such in- 
tensities have again a Poissonian distribution as the large time asymptotic, 
but with a constant intensity. This constant we can identify as the arithmetic 
mean of the initial intensity. For this result a careful consideration of the 
notion of convergence is required. This is one of the aforementioned techni- 
cal pitfalls caused by infinite volume. The above condition on the Fourier 
transform of the initial intensity cannot be directly described in terms of 
the intensity itself. It seems natural to claim that the large time asymptotic 
exists for all bounded intensities which have arithmetic mean. However, 
for example, bounded slowly oscillating intensities may not have arithmetic 
mean, cf. Remark 15.81 Under the condition of weak local asymptotic nor- 
mality of the one-particle process, in [DSS82J Theorem 4.1 the authors show 
that any free infinite particle dynamics is asymptotically similar to Brownian 
motion. We can show that that one-particle process is weak locally normal 
if the jump rate has finite second moments. Their proof is based on a cen- 
tral limit theorem type result and hence the existence of second moments is 
essential. The proof given in this paper is straight forward, simple and hold 
for general jump rates, but for a more restricted class of intensities. 

In Subsection 15.21 the so-called hydrodynamic limit for the free Kawasaki 
dynamics is considered. This means that the asymptotic of the corresponding 
scaled empirical field n[ e \if, X) := n(e d ip(e-), X e -« t ), k = 1,2, is the object 
under study. In order to obtain a nontrivial asymptotic, the initial intensity 
z has to be scaled as z(s-). The limiting process for e — > + is a deterministic 
evolving density profile given as the solution of a partial differential equation. 
For symmetric jump rates the obtained limiting equation, for n = 2, is a 
second order elliptic equation and the coefficients are given by the second 
moments of the jump distribution. For asymmetric jump rates, the limiting 



6 



differential equation, for k = 1, is of first order, cf. Proposition 15.151 and 
Remark 15.161 One may obtain a combined equation if one considers weak 
asymmetries, cf. Proposition 15 . 1 71 and Remark 15 .181 Actually, the limit gives 
a way to construct a solution for the partial differential equation. If the 
Fourier transform of the initial intensity is a measure, then these results can 
be obtained directly via Fourier transform, considering carefully the sense of 
convergence. If the jump rate has all moments finite, we are able to prove 
very strong convergence of the semi-group and hence can treat more general 
initial intensities. However, this is technically harder, because further careful 
approximations are necessary and a stronger sense of convergence has to be 
considered. The presentation of this result is postponed to Proposition 15.191 
for the clarity of the presentation. In [DSS82] the authors also consideres 
the Case 1 of Proposition 15.151 and 15.191 for jump rates with finite second 
moment. In Proposition 15.151 only finite first moments of the jump rate are 
required, but not all intensities are considered. 

The above described results about large time asymptotic and hydrody- 
namic limit are extended to general initial distributions in Section O We 
require some kind of mixing property for the initial distributions, the so- 
called decay of correlation, which we formulate in terms of the cumulants (or 
Ursell functions). For the hydrodynamic limit we have to require in addition 
that the first correlation function of the initial measure converges reasonably. 
The first correlation function is the density of the first moment measure of 
the empirical field. In the hydrodynamic limit the solutions of the associated 
partial differential equations are constructed as the limits of the first correla- 
tion functions, cf. Proposition 16. 11 The large time asymptotic is clearly again 
a Poisson measure for a constant intensity. However, the constant is now the 
arithmetic mean of the first correlation function of the initial distribution. 

We prove in Section I6T51 that the conditions required above are fulfilled, in 
particular, by Gibbs measures in the high temperature low activity regime. 

2 Kawasaki dynamics 
2.1 The generator 

The configuration space T := T R d over M. d , d G N, is defined as the set of all 
locally finite subsets of M. d , 

r := {7 C R d : |7a| < 00 for every compact A C R d ] , 
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where |-| denotes the cardinality of a set and 7a := 7flA. As usual we identify 
each 7 G T with the non-negative Radon measure ^2 xej 5 x G A4(R d ), where 
5 X is the Dirac measure with mass at x, ^ xg0 S x is, by definition, the zero 
measure, and Ai(M. d ) denotes the space of all non-negative Radon measures 
on the Borel a-algebra B(M. d ). This identification allows to endow T with the 
topology induced by the vague topology on Ai{R d ), i.e., the weakest topology 
on T with respect to which all mappings 

are continuous. Here C c (M, d ) denotes the set of all continuous functions on 
M. d with compact support. By B(T) we will denote the corresponding Borel 
a- algebra on T. 

Given a non negative function a G dx), the generator L := L a of 

the free Kawasaki dynamics for an infinite particle system is given by the 
informal expression 

(LF)( 7 ):=V/ dya{x-y)(F( 1 \xUy)-F(j)). (1) 

We proceed to give a rigorous meaning to the right-hand side of (0Q). Let 
O c (R d ) denote the set of all open sets in M. d with compact closure. A ^Im- 
measurable function F is called cylinder and exponentially bounded whenever 
there is a A G O c (R d ) such that F{i) = F( 7a ) for all 7 G T, and < 
Ce c ^, 7 G T, for some C, c > 0. For such a function F one has 



< 



J dya{x - y) |F( 7 \xUy)- F( 7 )| 





' dya(y) + \ 




m d J a 



a{x-y) 2Ce c(l7Al+1) , (2) 

which is finite provided the configuration 7 is an element in T a C T, 

r a := I 7 e r : y ^ ^ a{x - y) is Lj L oc (R d , dy) 

In this case, the sum and the integral in ([1]) are finite, and thus the operator 
L is well-defined on the space TL° eb (Y a ) of all cylinder functions exponentially 
bounded on Y restricted to Y a . 
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Concerning the set T a , we note that /x(r a ) = 1 for any probability measure 
H on T with first correlation function kjp bounded. This follows from the 
fact that for each closed ball B(n) C M. d centered at and radius n G N we 
have 

/ n(di) / dyS^a(x-y) = / dxkfi\x) / dya(x-y) 

Jr JB(n) xei JR d ' JB(n) 

< C\\a\\ L i {R d tdx) vo\(B(n)) 

for any constant C > \kft\. Here vo\(B(n)) denotes the volume of B(n) 
with respect to the Lebesgue measure on R d . Clearly, the latter implies that 

Mr\r a ) = o. 

In view of these considerations, throughout this work we shall restrict our 
setting to T a . 

Among the elements in jFLg b (r a ) we distinguish the functions e^(/), 
called Bogoliubov exponentials, 

e B (/,7):=n( 1 + ^ x ))' ^ GF ' 

for / being any bounded £>(lR d )-measurable function with compact support 
(/ G B c (M d )). A first reason to distinguish these functions is due to the 
especially simple form for the action of L on them, namely, for all / G B c (ML d ) 
and all 7 G T a , 

(Le s (/))( 7 ) = W dya(x-y)(f(y)-f(x))e B (f, 1 \x) 

xG 7 jRd 

= £(A/)(x)e B (/, 7 \x), (3) 

where 

Af(x):= I dya(x-y)(f(y)-f(x)). (4) 

In the sequel we call the linear operator A the one particle operator. Due 
to its special role throughout this work, its properties will be studied in more 
detail in the next subsection. 

Given a locally integrable function z > 0, we remind that the Poisson 
measure n z with intensity z is the unique probability measure on T for which 
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the Laplace transform is given by 



' vr,(d 7 ) e^' 7) = exp ( / dx {e^ x) - l) z{x) j , (5) 

for all ip in the Schwartz space V{M. d ) := C£°(R d ) of all infinitely differen- 
tiable functions with compact support. We recall that a measure \i on T a is 
called infinitesimally reversible with respect to the operator L defined in ([I]) , 
whenever L is symmetric in L 2 (T a ,fi). 

Lemma 2.1 Assume that a is an even function. Then for any real number 
z > the Poisson measure tc z with constant intensity z is an infinitesimally 
reversible measure with respect to L. 

Proof. For all F, G G FL%(r a ) we have 

/ ^(^7)^(7)^(7) = / 7r,(rf 7 )V/ dya{y-x)F{ 1 \x\Jy)G{ 1 ) 

Kz(d-f)Y] / dya{x-y)F{j)G{j). 



The result follows by applying twice the Mecke identity (cf. |Mec67t Theo- 
rem 3.1]) to the first integral. ■ 

Remark 2.2 It is clear that the linear space spanned by the class of func- 
tions es(/), / G B c (R d ), is in jFLg b (r a ). The space jFL° fe (r a ) also contains 
the class of coherent states e^Xf) corresponding to / G B c 



e n z (f) ■= ex P y~ J ^dxz(x)f(x)j e B (f). (6) 
2.2 The one particle operator 

In the sequel the one particle operator A introduced in (TJJ will play an 
essential role. Because of this, in this subsection we shall collect its main 
properties used below. We observe that in stochastic analysis the operator A 
is known as a relatively simple example of a bounded generator of a Markov 
jump process on M. d (see e.g. |EK86[ Section 4.2]). 
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In terms of the usual convolution * of functions, we also note that 



Af = a*f- a {0) f, a {0) := / a(x) dx. (7) 

Therefore, the properties of convolution of functions (namely, Young's in- 
equality) lead straightforwardly to the L p results stated in the next propo- 
sition. There, we also consider the real Banach spaces B(M. d ) and C 00 (]R rf ), 
respectively, of all bounded measurable functions and of all continuous func- 
tions vanishing at infinity, both with the supremum norm ||/|| u := sup^^d \f(x) 
We recall that a strongly continuous contraction semigroup (T t ) t > is called 
sub-Markovian whenever for < / < 1 it follows < T t f < 1 for every t > 0. 
If, in addition, T t f n /* 1 for some sequence /„ / 1, then (T t ) t > is called 
Markovian. Here, for the spaces B(M. d ) and C 00 (M d ) the convergence is point- 
wise, and for an L p -space the convergence is almost everywhere. A strongly 
continuous contraction semigroup (T t ) t > is called positivity preserving, if 
/ > implies T t f > for every t > 0. 

Proposition 2.3 The linear operator A is a bounded operator on L p (M. d , zdx), 
for z > constant and 1 < p < oo (on B(M. d ) and on Co^lB^) ). As a conse- 
quence, A is the generator of a uniformly continuous semigroup (e tA ) t >o on 
L p (R d , zdx), 1 < p < oo (on B(R d ) and on Coo(M d ) ), 



e tA. ( tA ) n 

n=0 



the sum converging in norm for every t > 0. Moreover, (e tA )t>o is a positivity 
preserving (contraction) semigroup on each L p (M. d , zdx) space, 1 < p < oo 
(on B(R. d ) and on C^W 1 ) ), Markovian on L p (IR d , zdx) for 1 < p < oo 

For the proof see e.g. [EK86[ Section 4.2], |Jac01| . 

Due to (J7j), the operator A, as well as the semigroup (e tA ) t >o, both either 
on a L p (M. d , zdx) space, 1 < p < oo, or on C^IR^), may be expressed in 
terms of the Fourier transform, 



see e.g. |Jac01j . 



11 



Proposition 2.4 For every function if in the Schwartz space S(R d ) of tem- 
pered test functions one has 

Ap{k) = (2vr) d/2 (o(Jfe) - a(0)) <f(k), k G R d , 

and 

! / .-/I. ™\ <Y"W\d/2/ 



(e M ^) (x) = -rA-nz / dke^e*^ 1 ^~^0(k), xeR d . (9) 
(27r) d / 2 J Rd 

Remark 2.5 Since a is non-negative, for all k G R d one has 

Re(a(Jfe) - a(0)) = 1 / rfx (cos((Jfe, x)) - l)a(z) < 0. (10) 



(270^,^ 

The equality in ( flOi) holds only for = 0. This follows from the fact that the 
set {x : (k, x) — 2mr, n G Z} has zero Lebesgue measure if and only if k ^ 0. 

Remark 2.6 According to Proposition I2.4[ the semigroup (e tA ) t >o is defined 
by a kernel /i t , t > 0, whose Fourier transform is given by 



3 t(27r) d / 2 (a(fc)-a(0)) 



(2tt)<V 2 

Because a is non-negative, a and thus /t f , t > 0, are positive definite functions. 
Through Bochner's theorem, the latter means that each /i t is a non-negative 
finite measure on R d . We note that the Markovian property of {e tA ) t >o 
means that each \x t is actually a probability measure. For more details see 
e.g. |Jac01j . 

Remark 2.7 It is easy to check that on the space L 2 (R d , zdx) the adjoint 
operator A* of A is defined by 



(A*f)(x) = / dya(y-x)(f(y)-f(x)) + / dy(a(y)-a(-y)) f(x), f G L 2 (R a , zdx). 

As a consequence, if a is an even function, then A is a bounded self-adjoint 
operator on L 2 (R d , zdx). In this case, it follows from (|SJ) that the semigroup 
(e M )t>o is a self-adjoint contraction on L 2 (R d , zdx). 

Remark 2.8 For a non-constant activity parameter z, the operator A is 
still bounded but, in general, the semigroup (e tA ) t >o is not any longer a 
contraction. For instance, on R consider a(x) = e~ x and z(x) = 1 + 
e -x +4a;_ ^ simple calculation shows that for <^(x) = 7r^ 1 / 2 e _z one has 
L da; <f(x)(A(p)(x)z(x) > 0, proving that the semigroup (e tA ) t >o cannot be a 
contraction in L 2 (R d , zdx). 
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3 Independent infinite particle processes 



3.1 A probabilistic approach 

As we mentioned at the beginning of Subsection I2.2[ the operator A is the 
generator of a Markov jump process on M. d . Following e.g. [KK86], we can ex- 
plicitly construct this process as follows. Consider the Markov chain (y^eNo 
on R d , with transition density function fi(x,y) = ■ Let (Z t )t>o be a 

Poisson process with parameter independent of the Markov chain. We 
then define the Markov process (X t )t>o by X t := Yz t , t > 0. This process 
has generator A and, by construction, it has cadlag paths in ~R d . We denote 
by D([0, oo), M. d ) the set of all cadlag paths from [0, oo) to M d and by P x the 
path-space measure corresponding to the process X starting at x G M. d . By 
E x we denote the expectation w.r.t. this measure. 

The process on V corresponding to L is the following random evolution: 
each particle evolves according to the above jump process, independently of 
the other particles, cf. Lemma 13.11 below. This independent infinite particle 
process was rigorously constructed in |KLR09j (see also [KLR03] ). In the 
following we present the main results therein. 

Assume that there exist an a > 2d and a C > such that 

C 

< a(y) < t 1— rr— , for all y G R d . (11) 

The condition sufficient for the construction done in [KLR09j is then fulfilled. 
Hence there exists a Markov process (-D([0, oo), 6), (X. t )t>o, (P-y)-yee) on the 
set of all 7 G T such that, for some m G N (depending on 7), 

|7B(n)| < mvol(S(n)), Wn G N. 

Here D([0, 00), 0) denotes the set of all cadlag paths from [0, 00) to 0, the 
process X t : D([0, 00), 0) — > is the canonical one, i.e., X t (u;) := w(t), 
uj G D([0, 00), 0), and each P 7 is the path-space measure of the process 
starting at a 7 G 0. By E 7 we denote the expectation w.r.t. P 7 . 

Note that the process cannot start at any arbitrary initial configura- 
tion 7 G T, as follows implicitly from the discussion before Theorem 2.2 in 
[KLR09J. One is obliged to restrict the set of possible initial configurations 
e.g. to 0. However, by |KKK04j . [Kun 99j. one has /x(0) = 1 for every prob- 
ability measure fi on T whose correlation functions k^, n G N, fulfill the 
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so-called Ruelle bound, i.e., there is a C > such that < C n for every 
n G N. This holds, for instance, for Gibbs measures w.r.t. superstable, lower 
and upper regular potentials, cf. |Ruc70j. In particular, it holds for Pois- 
son measures with intensity being any non-negative bounded measurable 
function. For a constant intensity this result was shown using ergodicity in 
[NZ76] . 

Choosing, for a fixed initial configuration 7 G 6, an enumeration {x n } ng N, 
the infinite particle process can be described more explicitly. For each 
n let us consider an independent copy of the one-particle jump process 
((X^) t > , P Xn ) introduced at the beginning of the section. In [KLR09] it 
is shown that a.s. the sequence (X^ ) n eN has no accumulation points and 
no two points in it which coincide. Then ({X^} ne ^) t > is the correspond- 
ing process on the configuration space 6 and the path-space measure is the 
"symmetrization" of (£)™ =1 P Xn - Moreover, the transition probability (Pt)t>o 
of the process (K t )t>o is just the product of the one-particle transition prob- 
abilities e tA (x — y)dy, i.e., YlnLi etA ( x n — y n )dy n - As a consequence, for all 
non-positive (p G V(M. d ) we have 

E 7 [e^] = j P t ( 7 , dO e<^> = J[E X [e**>] , (12) 

cf. |KLR09j (see also Lemma 1X21 below) . or in terms of Bogoliubov exponen- 
tials for ip G V(R d ) with -1 < if < 0, 

E 7 [e B (<f,X t )} = \{E X [(p(X t ) + 1] = e B (e^, 7 ). (13) 

£€7 

More generally, it follows that 

E 7 [ e ^^ • . . e ^- x O] =J[E x [e^ x ^ ■ ■ -e^ x ^] , 

for all < t\ < . . . < t n and all non-positive (f\, . . . , (p n G T>(Mr), n > 2. By 
a monotone approximation procedure using Riemannian sums, this relation 
allows us to calculate the Laplace transform of the path-space measure P 7 , 
i.e., 



E 7 



-/dt<VJ(t,-),X t ) 



l[E x [e-^ dt ^ Xt ^, (14) 



XG7 
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for all non-negative continuous functions ip from [0, oo) x M. d to K with com- 
pact support. 

The next result gives a relation between the process X and L. 



Lemma 3.1 For all F G jFLg 6 (0) and all 7 G 6 it holds 



F(X t ) - F(X ) - / dsLF(X s 
Jo 



0. 



(15) 



Proof. For all F of the form esif) with ip G D(IR d ) one has 



^E 7 [F(X«)] = 4 B (e''V, 7 ) =£* u # t (e'V7\4 



XG7 



Hence using the product structure of the path-space measure and 
finds 



one 



^Ae tA p{x)e B {e tA p^\ 



x£7 



Av?(a;)e i3 (v3,Xi\ x) 



E 7 [Le B (^)(X t 



In order to make the previous calculations rigorous and to extend them 
to the whole space JFLg b (6), we have to establish bounds for the considered 
expressions. According to for all F G JFL° 6 (6), one may bound the 
integrand LF(X S ) in (fT5l) by 



|x s nA| j 


dya(y) + / 




R d J A 



a(x - y) 



xex s 



2C < eC (|x s nA|+i) > 



For some constants C, C" > one may bound the previous expression by 
C"e^ x ' Xs \ where x(v) = n^pifiwa with a as in ffTTj) . The first summand in ffT5l) 
and the expression XLex t I^V 9 ( a; )l e B(lv 9 l) can be bounded analogously. 

Due to Lemma 13.21 below the expectation of these bounds are finite for all 
7 G 9. ■ 



c 



(1+12/ 



lyjya, for some C > 



and a > as in <f771). 



Lemma 3.2 Lei 

T/ien £/iere exists a c > sitc/i i/ia£ A% < c% and e x < eC *X / or all t > 0. 
Moreover, for all 7 G B and a// measurable functions <p such that ^ is 
bounded, one has that the product f^^g (1 + V 9 ( x )) absolutely convergent 
and 

E 7 [e (v ' Xf> ] < 00. 
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Proof. The previous considerations yield 

E 7 p' Xt >] < E 7 [A Xt >] = l[E x [e^] = ]Je tA e^x). 

2G7 

Thus the proof amounts to show the convergence of the latter infinite product. 
For this purpose, we note that 

^ x " s) = (i + 1, 5 a $+\ff = (1 + Wr/2 * W ' 

This gives for the one particle operator 

A X {x) = [ dya(y)x(x-y)-a (0) x(x) < ( [ dy{l+\y\) a ' 2 a{y))x{x)-a^x{x). 
Jm d v JB, d J 

Due to the bound (iTTj) for a, the integral J Rd dy(l + |2/|) a ^ 2 a(|/) is finite. 
By Gronwall's lemma this yields E x [x(X t )} = e tA x(x) < x(x)e rf , where c = 
/ Rd dy(l + \y\) a/2 a(y)-a(°\ Therefore, also |e' A (e x -l)| decays as (l + \y\)~ a/2 . 
According to the definition of G, one can prove that (x, 7) < 00 for all 7 G 0. 
Hence also the product Ylxej e tA e x {x) converges absolutely and is finite. ■ 

In Sections and one considers the one-dimensional distributions of 
processes starting with initial distributions /i which are probability measures 
on 0. The path-space measure P M corresponding to such a process is given 
by J* e P 1 ii(d'y). Its one-dimensional distribution is a probability measure P* t 
on defined for all non-negative measurable functions F by 

/ P*{dl)F{l) := / M^7)E 7 [F(X t )] . (16) 
Je Je 

In particular, for /x = S 7 , 7 G 0, the one- dimensional distribution coincides 
with the transition kernel Pt(7, •) described above. 

For functions F being Bogoliubov exponentials, definition (|T6j) leads to 
the so-called Bogoliubov functionals |Bog46| . By definition, the Bogoliubov 
functional corresponding to a probability measure /1 on T is defined by 

:= J e B &,lMdj), V e V(R d ). 

Such a functional is an analogue of the Fourier-Laplace transform on configu- 
ration spaces, cf . |Kun09j . Due to (fT5|) there is an interesting relation between 
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the Bogoliubov functional corresponding to the initial distribution and the 
Bogoliubov functional corresponding to the one- dimensional distribution of 
the process at a time t > 0, namely, 

P*(d7)es(y?,7) = / Kdl)e B (e tA <p n ) = B,(e tA cp). (17) 
e Je 

In particular, for /i = tt z for some bounded intensity function z > one 
finds 

e B (^7)P* t (d 7 )= / e B (e tA <p,i)TC,(d'y) = ex P ( f e tA V {x)z{x)dx\l$) 
'e Je \JR d J 

for all if G X>(R ). In Section [5] we shall consider this special case in more 
detail. 



3.2 An analytic approach 

Within the framework of infinite dimensional analysis on configuration spaces 
|KK02j one may derive alternative representations and constructions of the 
dynamics. Instead of describing the infinite particle dynamics on T through 
the Kolmogorov equation -^F t = LF t , the so-called if-transform |Len75at 
ILen75bj allows an alternative description for the action of L on TL° eb {T a ). 

Given the space r := {7 G V : |7| < 00} of finite configurations en- 
dowed with the metrizable topology as described in |KK02] , let B exp j s (T ) 
be the space of all exponentially bounded Borel measurable functions G (i.e., 
|G(^)| < C\e° 2 ^ for some Ci,C 2 > 0) with local support (i.e., there is 
a A G O c {R d ) such that G ta e ro:|»jn« d \A)|^o}= 0)- The f^-transform of a 
G G B exp j s (T ) is the mapping KG : T — > R defined for all 7 G T by 
{KG){i) := £ ^ Gfo). Note that K (B exp , ls (T )) C Given a 

|tj| <oo 

G G -BezpJstTo)} then in terms of the operator L we obtain 

{L(KG))(i) = J2 f dya(x-y)[(KG)( 1 \x) + (KG(-Uy))( 1 \x) 

xei J»*\<a\*) 

-{KG){ 1 \x)-{KG{-Ux)){ 1 \x)\ 
which leads to the so-called symbol acting on quasi-observables L, 

(LG)(rj) :=V / dya(x - y) (G(r, \ x U y) - G{rj)) , G G iW*( r o), 
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and the corresponding time evolution equation -j^G t = LG t - The transition 
kernel P t corresponding to L is then given by 

I P t (r)',dr))G(r))= I f[ [ dy t e tA {x % - y t )G({ yi , . . . , y n }). 
Jr JK. dn i=1 Jm. d 

This allows us to extend the explicit formula for transition kernels of X to 
the class of all so-called observables of additive type 

jf P t ( 7 , dfl (KG)(0 = K Qf P t (; drj) G(r,)j (7). 

By duality one may extend the dynamical description to correlation func- 
tions. For this purpose, on the space Tq = U^L {7 G T : \^y\ = n} let us con- 
sider the so-called Lebesgue-Poisson measure \ z := Yl^Lo "1 { zm )^ a \ where 
m denotes the Lebesgue measure and each (zm)^, n G N, is the image 
measure on {7 G Y : I7I = n} of the product measure z{x\)dx\ ■ • ■ z(x n )dx n 
under the mapping (M, d ) n 3 (xi,...,x n ) t— > {xi, x n }. If for some proba- 
bility measure fi on T there is a function on T such that the equality 
f r fi(dj) (KG) (7) = / ro \(dr)) G(r])k^(r]) holds for all G G B expM (Y Q ) we call 

the correlation function corresponding to fi. Here we abbreviate A := Ai. 
Denoting by L* the dual operator of L in the sense 

/ X z (d V ) (LG)( V )k( V ) = [ X t (dr,)G(r,)(L*k)(r,), 

one obtains the following expression 

„.<-„ JR d 



The corresponding time evolution equation for correlation functions, ^k t — 
L*k t , is the analogue of the BBGKY-hierarchy for the case of the free Kawasaki 
dynamics. In our case this equation can be explicitly solved, namely, 

» n 

kt({xi,---,x n })= dy 1 ---dy n k f ,({y 1 ,---,y n })Y[e tA (y i -Xi) (19) 

for the initial condition k^. Let us note that if one assumes a Ruelle bound 
for the initial correlation function k^, then all the above considerations can 
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be made rigorous. Moreover, similar arguments used in [KKZ06J show that 
each k t is actually a correlation measure corresponding to some probability 
measure fit on T. 

According to the previous considerations, the time evolution of the parti- 
cle system may also be described in terms of Bogoliubov functionals, cf. [KKO 06J. 
through the time evolution equation 

7zBM=f dx l dya(x-y)(^y)-<p(x)) 5 -^^-, <p G V(R. d ). 
at J R d j Rd 6cp{x) 

Here SB S ^'^ denotes the first variational derivative of B fl:t at <p. Actually, 
Bnt((p) is the Bogoliubov functional corresponding to the one-dimensional 
distribution of the process starting in fi. Hence, the previous equation has 
an explicit solution given by (fTTl) . i.e., B^f) = ^ti{ e f)- 

4 Equilibrium dynamics 

In this section we are interested in the representation of the generator L 
and its semigroup in terms of the creation, annihilation and second quanti- 
zation operators as analytic expressions. These representations are possible 
because there is a well-known canonical unitary isomorphism between the 
(symmetric) Fock space and the Poisson space. We start by recalling this 
isomorphism. Our approach is based on [AKR98] and [KSSU98] . but see also 
[KKO02J and references therein. 

4.1 The L 2 (tt z ) space and the Fock space representation 

Let z be a non- negative constant. We consider the complex Hilbert space 
L 2 (n z ) := L 2 (T, B(T), n z ) of square integrable complex valued functions on 
T with respect to the Poisson measure ir z . The coherent states introduced in 
Remark 12.21 generate the system of so-called Charlier polynomials, namely, 

n=0 

where (T> f )^ n is the n-th symmetric tensor product of the Schwartz distribu- 
tions space V'{R d ). This system is orthogonal and any F G L 2 (n z ) can be 
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expanded in terms of Charlier polynomials 

oo 

Hi) = j>£(7), f (n) ^ f (n) e L \ z dxt n . (20) 

n=0 

This yields a unitary isomorphism I 7Tz between L 2 (tc z ) and the so-called 
symmetric Fock space 

oo 

T(L 2 {zdx)) := n\L 2 {zdxf n , L 2 {zdxf° := C. 

n=0 

More precisely, for each F 6 L 2 (n z ) of the form (120]) one has I Vz {F) = 
(f^)^Lo- Next we recall the definition of annihilation, creation and second 
quantization operators on the total subset of Fock space vectors of the form 
j{n) _ _ _ _ <g>/ n . The action of the annihilation operator a~(h) of a h E 
L 2 (zdx) on f( n ' is given by 

n 

a-(h)fW : = J2(h, ■ ■ ■ ■ ■ ■ ®/„ G L^xf^. (21) 

3=1 

The adjoint operator of a~(h), called creation operator and denoted by a + (h), 
acts on elements /< n ) e L 2 {zdx)® n by a + (h)f^ = h®f n \ 

Given a contraction semigroup (e tA ) t >o on L 2 (zdx) one can construct a 
contraction semigroup (Exp(e* A ))t>o on F(L 2 (zdx)) defined by e tA (8) ... ® e M 
on each space L 2 (zdx)® n . Its generator is the so-called second quantization 
operator dExpA corresponding to A. Hence the image of the Fock coherent 
state e(f) := U'® n /n\)™ =Q , f G L 2 (zdx), under Exp(e M ) is given by 

Exp(e M )(e(/)) = e(e M /). (22) 

Through the unitary isomorphism I nz we obtain a contraction semigroup 
(Exp 7rz (e M )) t >o on L 2 (n z ). In particular, since I~ 1 e(f) = e Wz (f), it follows 
from (jf)]) that 

Exp^(e M )e B (/)=e B (e* A /). (23) 

In our case, for non-constant functions z the semi-group e is in general not 
a contraction, see Remark 12.81 and indeed the previous construction ought to 
fail, cf. Remark 15. II 

Finally, we would like to present the "annihilation and creation operators" 
in the form more common in physical literature. In this heuristic way one can 
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also treat non-constant intensities. For each x G M. d we define an operator 
a~(x) acting on / = (/ (n) )„ by 

(a-(x)f)W(y u ...,y n ) = V^+Tf n+1 \x, yi , . . .,y n ). 

The adjoint of the operator a~(x) is formally given by 

1 n 

where y& means that the fc-th coordinate is excluded. Actually, the expression 
for a + (x) is well-defined as a quadratic form. It is easy to check the relations 
a ± (f) = J Rd a ± (x)f(x)dx in the sense of quadratic forms. 

The operator L can be expressed in terms of creation and annihilation 
operators, namely, 

L = dEx Pwz (A) + a- z (A* z l), (24) 

where A* is the adjoint operator of A with respect to the scalar product on 
L 2 (R d ,z{x)dx), dExp (A) and a n are, respectively, the image of dExp(A) 
and a~ under 7^. For more details see e.g. [KKO0 2j|. Note that the expres- 
sion of the annihilation operator a~ z in L 2 (7r 2 ) depends on the intensity of 
the underlying Poisson measure. Rewritten in a style more similar to the 
common usage in physics, (1241 takes the form 

J dxz(x) J dy (a(x - y) - a (0) S(x - y)) (a+(x)a~(y) - a~(y)) , 

where a^(x) are the image of a ± (x) under 7^. 

4.2 The symmetric case 

If a is an even function, the situation is essentially simpler (see Lemma 12.11) 
and yields an alternative construction to the one presented in Section [3] As 
a matter of fact, for a an even function, and for each constant z > 0, the 
operator L defined in (DO) gives rise to a Dirichlet form on L 2 (T, tt z ), 

[ v z {d 1 ){LF){ 1 )F{ 1 ) = -\ /vr,(rf 7 )V / dya(x-y) \F( 7 \x U y) - F( 7 )| 2 . 
Jr 1 Jr xey^ Rd 
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This allows the use of Dirichlet forms techniques to derive a Markov process 
on r with cadlag paths and having ir z as an invariant measure |KLR07j . Ac- 
tually, one can show that in this situation L is the second quantization opera- 
tor corresponding to the non-positive self-adjoint operator A on L 2 (M. d , zdx) 
(Remark 12.71) . Hence, L is a negative essentially self-adjoint operator on 
L 2 (T, 7r 2 ), and it is the generator of a contraction semigroup on L 2 (T, tc z ). 

4.3 Second quantization operators 

According to (fl3l) and ([6]) the action of the transition kernel on coherent 
states corresponding tt z is given by 

E 7 [e nz (<P,Xt)} = e^(e M ^, 7 )exp (j^dx (e tA ip{x) - <p(x))z\ . (25) 

This allows to express the action of the transition probability on coherent 
states in terms of annihilation and creation operators 

Exp^(e iA )e a - (A * 1) . (26) 

Observe that the right-hand side of (j2"51) gives the action of a semigroup 
which preserves coherent states. Lemma 14.11 shows that L is the generator 
of a strongly continuous Markov semigroup. 

Lemma 4.1 Let P co h be the vector space spanned by all functions esi'-p) with 
ip G L 1 (M d , zdx) fl L 2 (R d , zdx). The operator L restricted to "D co h is closable 
in L 2 (T,it z ) and its closure is an extension of the operator (L, !FL^ h (T a )) 
defined in Section [H Moreover, it is the generator of a Markov semigroup 
and e tL = Exp nz (e tA ). 

The proof is an adaptation for a non-symmetric generator of the proof 
of Proposition 4.1 in [AKR98J. Technically, in our case the proof is simpler, 
because A is a bounded operator and we consider functions of the type es{(p)- 

Proof. According to the bound (T5]), one can calculate the action of the 
adjoint of (L, jFLg b (r a )) by the Mecke formula, i.e., 

(L*F)( 7 ) := Y, I d V («(*/ - X ) F ^ \xUy)-a{x- y)F(>y)) , VF G ^L c ° b (r a ). 
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Hence L* is densely denned and L is closable on any dense subset of J r L^, h (T a ). 
Using again the bound one sees that the closures of (L, T> coh ) and (L, jFL° b (r a )) 
coincide. As A is a bounded operator in L 1 (M d , zdx) and in L 2 (M. d , zdx), the 
space T>(M. d ) is a core of A in L 2 as well as in L 1 . Define an L 2 (T,tt z ) 
semigroup on £> co h (as an extension by continuity) of the following explicit 
formula 

T t e B (<p) := e B (e tA <p). 

Hence £> C oh is a core for the generator of T t . By computation one sees that 
if i — > gb{(p) is a differentiable function from L l (M. d ,zdx) n L 2 (M d , zdx) into 
L 2 (r,7r z ) with derivative given by 

Summarizing, this yields that T t is a strongly continuous contraction semi- 
group on L 2 (r, 7i z ) such that for all <p G L 1 ^, zdx) n L 2 (M d , zrfx) holds 

^T t eBfoO(7) = E ^^(x)e B (e 4A ¥,, 7 \ x). 

For </? G the r.h.s. is just Les(e tA (p). Since A is bounded this can be 

extended to all ip G L 1 (M. d ,zdx) H L 2 (R d , zdx). These two facts combined 
yield e tL = T t and V coh is a core of L. Moreover, the operator L coincides 
with the second quantization operator dExp Wz (A) on V coh . This space is also 
a core for dExp nz (A), cf. [BK88]. Hence the two operators coincide. ■ 

5 Local equilibrium dynamics 

The previous considerations were mainly for Poisson measures 7r 2 which have 
as an activity parameter a constant function z. According to Lemma 12.11 
and Section [21 such measures are reversible measures for the free Kawasaki 
process. As a first step towards other initial distributions, in this section we 
consider the so-called local equilibrium case, that are, Poisson measures with 
non constant activity parameter. In this case, if the activity z is a slowly 
varying function, then in a small volume the Poisson measure effectively has 
a constant activity. This property justifies the chosen name for this section. 
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Remark 5.1 Although for bounded non-constant functions z the process 
starting in n z can be constructed, cf. Subsection 13.11 one cannot expect that 
the expression given in (fTTj) can be extended to a semigroup either in the 
L l {T a , 7t z ) or in the L 2 (T a ,7r z ) sense. The existence of the semigroup in an 
L p -sense can only be expected w.r.t. an invariant measure. This can be seen 
from the following equality 

e tL e B (f)\\„ M , . , ,,, M , ,..,,„ ,, „ z{x)dx 



\\e B U)\\Lv^ z ) 
and the fact that 



expfi / (\l + e tA f(x)\ p -\l + f(x)\A, 
\P Jr* v j 



\\e tL \\ 0vLP(7T) > sup expf-/ f(l + e tA f(x)r-(l + f(x)r)z(x)dx 

^' f£V (R d ),/>0, \P JR d v 7 

where the r.h.s. is infinite if (e tA )t>o is not a contraction semigroup (see 
Remark El]). 

Let z > be a bounded measurable function. According to (1181) . for the 
one-dimensional distribution of the free Kawasaki process (X t ) t >o with initial 
distribution tc z one finds 

e B ^,l)P*t(dj)=exp( [ dx e tA <p(x)z(x)) , (27) 

for every ip £ T>(M. d ) and every t > 0. 

For each t > fixed, let us now consider the linear functional defined on 
L 1 (M <i , dx) by 

L\R d , dx) 3 f h-> / dx(e iA /)(x),s(x). 

</R d 

Due to the contractivity property of the semigroup e* A in L 1 (M d , dx), see 
Proposition 12.31 and to the boundedness of z, this functional is bounded on 
L 1 (M d , dx), and thus it is defined by a kernel z t £ L 00 (R d , dx), that is, 

dxe tA f(x)z(x) = / dxf(x)z t (x), (28) 

JR d 

for all / £ L 1 (M d , dx). Moreover, since e tA is positivity preserving in L 1 (IR c( , dx), 
it follows from ([28]) that z t > 0. 
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This shows by ([27]) that the one- dimensional distribution P* t is just the 
Poisson measure n Zt (see also [Dob56j . |Doo53j ). 

Furthermore, the path-space measure P 7 ^ corresponding to the process 
is also Poissonian. It is a Poisson measure on Pdqo.oo),^) with intensity P z , 
where P z is the path-space measure on D([0, oo), M. d ) of the one-particle jump 
process corresponding to A with initial distribution z(x)dx. Using the fact 
that the path-space measure P 7 is supported on D([0, oo), 6) one easily sees 
that this implies that P nz is actually supported on a subset of r_D([o,oo),R d ) 
which can be naturally identified with D([0, oo), 0). 

Lemma 5.2 The path-space measure P nz of the process X starting with ini- 
tial distribution tc z is the Poisson measure npz on T D(\Q,oo),R d )- ^ n particular, 
for any continuous function if on [0, oo) x M. d with compact support we have 



P*'(dw)e-/o"M*.0M*)>«« = exp ( / e -/ ~ *>(*,«(*))* _ i 

Proof. By the definition of the path-space measure corresponding to the 
initial measure 7r z one has 



V^{du:)e 



D([0,oo),: 



e 



TT z (d^). 



Due to pijl the latter expectation is equal to 



e B (E. 



e 



-f™ v (t,X t )dt 



-l,j)-K z (d-f) = exp 



E , 



z(x)dx 



which yields the required result. 



5.1 Large time asymptotic 

In this subsection we want to study the behavior of the one-dimensional 
distribution of the free Kawasaki dynamics for large times. As mentioned 
above (1281) . the free Kawasaki dynamic leaves the Poissonian structure of 
the initial distribution unchanged and only the underlying intensity develops 
with time. Thus, the analysis of the large time asymptotic behavior reduces 
to the large time asymptotic of the intensity, cf . Lemma 15.91 In other words, 
we reduced the problem to a one particle problem where the intensity plays 
the role of the initial distribution. In the following remark we discuss which 
initial distributions are natural in this context. 
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Remark 5.3 As we work in the framework of a system scale much larger 
than the time scale, namely, we first perform the thermodynamic limit and 
only afterwards consider the time asymptotic, the situation is more subtle 
than it might seem at a first glance. The question is which class of initial 
distributions is the natural one. Usually, in the study of an one particle 
system one assumes that the initial data is integrable. In this case the dy- 
namics is ergodic and the probability measure concentrated on the empty 
configuration is the invariant measure. Physically, this situation describes 
systems with zero density and perturbations of them. These perturbations 
are already singular in the sense that the corresponding Poisson measures 
are mutually singular. From a physical point of view, non-zero densities are 
interesting. The corresponding invariant measures are the Poisson measures 
with constant intensities z > 0. However, these intensities are not any longer 
integrable perturbations of the zero intensity. Furthermore, also their mu- 
tual differences are not integrable. Hence a natural setting for the initial 
intensities seems to be bounded non-negative measurable functions. 

Definition 5.4 One says that a function z 6 Lj 0C (R d , dx) has arithmetic 
mean, denoted by mean (z), whenever the following limit exists 



Note that not every bounded function has arithmetic mean, cf. Remark 1 5. 8 1 
The following proposition states that the large time asymptotic of the one 

particle distribution is still Poissonian with intensity given by the arithmetic 

mean of the initial intensity. 

Proposition 5.5 Let z > be a bounded measurable function whose Fourier 
transform is a signed measure. Then z has arithmetic mean and the one- 
dimensional distribution t converges weakly to 7r mean(z ) when t goes to 



This result is proved combining Lemma 15.91 and Corollary 15.121 
Note that, in particular, one sees that the process is not ergodic, because 
the large time asymptotic depends on the initial distribution. The situation 
of integrable intensities is technically comparable with the case in which the 
time scale is much larger than the space scale, e.g., instead of M. d one works 
on a torus. 




(29) 



infinity. 
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Under the assumptions of Proposition 15.51 the activity z has arithmetic 
mean, cf. Corollary 15. 121 Then, using Fourier transform, one can easily derive 
the large time asymptotic, cf. Lemma [5.101 The same technique also yields 
that the large time asymptotic only depends on the space asymptotic of the 
intensity z, cf. Corollary 15.141 

In Proposition I5.5[ the assumption concerning the Fourier transform is 
actually not elegant, because we cannot reasonably restate it in terms of z in 
the position variables. Therefore, we give several examples below, cf. Exam- 
ple 15.71 and we derive certain properties of the arithmetic mean. Some types 
of reasonable asymptotic behavior do not fulfill the aforementioned Fourier 
transform assumption, e.g. Example I5.7( v). 

Remark 5.6 If two functions Zi,z% G L 1 1 oc (IR d , dx) have arithmetic mean, 
then for every ai,«2 6l the function a\Z\ + a^z^ also has arithmetic mean 
and mean(ai,2i + 0:2-22) = O'imean(zi) + a^mean^). 

Example 5.7 To be more concrete we give some examples: 

(i) If z is a constant function then mean(2;) = z. 

(ii) If z decays to zero, i.e., for every e > there exists an R > such that 
\z(x)\ < e for x £ B(R), then mean(z) = 0. 

(iii) If z G L p (R d , dx), p G [1, 00), then mean(z) = 0. 

(iv) Also trigonometric functions have no influence, e.g., for d — 1 and 
z(x) = 1 + esin(x) we have mean(2;) = 1. 

(v) A less trivial example is the following one. Given z ,zi > 0, let z be 
the function defined at each x — (xi, . . . , x<i) G M. d by 

/ x _ J Z\ if X\ > 
\ z otherwise 

In this case one finds mean (z) = Z °+ Zl . Note that in this case z is not 
a signed measure. 

Remark 5.8 The arithmetic mean does not exist for all bounded non-negative 
functions. 
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1. For d = 1 and 



z[x) 



1 if 2 2k < \x\ < 2 2k+1 for a k G N 
otherwise 



the integral in (1291) oscillates between 1/3 and 2/3. Thus z does not 
have arithmetic mean. 

2. Another example is given by the following slowly oscillating function 
z(x) = cos(ln(l + \x\)) + z Q , x G R d , 
where is a constant greater or equal to 1. Then for large R it holds 

i/p/pNN / rfx ~ 7j== to sin (M R + 1 ) + arctan(d) ) + 2r . 

vo\{B(R)) J b{R) Vl + d 2 V / 

In general slowly varying functions will show in general spurious be- 
havior. 

First, we prove that one may reduce the proof of Proposition 15.51 to a one 
particle system. This follows from the independent movement of the particles 
discussed in Section [3j 

Lemma 5.9 Let < z G L 1 1 oc (M d , dx) be such that 

lim / dxe tA (f(x) z(x) =: / z 00 (dx)(p(x) 

t^+coJ Rd J Rd 

exists for all <p G V(R d ). If (non-negative) Radon measure on ~R d , 

then the one- dimensional distribution t converges weakly to tc Zoo when t 
tends to infinity. 

Proof. According to (J27J), for all <p G V(K d ) we have 

/ p *,t( d ^ e(v ' 7> = ex P (/ dxe u (e v - l)(x) z(x) 

with — 1 G X?(IR d ) too. By assumption, the latter converges when t — ► oo 
to 



cxp 



Zoo (dx) (e^ - 1) 
28 



By the definition (jSJ) of a Poisson measure, observe that this is the Laplace 
transform of ir Zoo . The convergence of Laplace transforms implies weak con- 
vergence, see |Kal76t Theorem 4.2]). ■ 
It remains to derive the large time asymptotic for the one particle system. 

Lemma 5.10 Let z > be a bounded measurable function such that its 
Fourier transform z is a signed measure. For all <p G 5(M d ) one has 

lim / dxe tA p{x) z{x) = (2tt)^ /2 z({0}) / dxp{x). (30) 

t^+oo J Rd - J Rd 

Proof. Through the Parseval formula and the explicit formula ([9]) for the 
semigroup (e tA ) t >o one finds 

dxe tA if(x)z(x) = [ z(dk)e t{2n)d/2 ^- k) - &m <p(-k). 



Since for all k G R d , Re(a(k) - a(0)) < (cf. Remark [23]), for every t > 
we have 

e t{2K)*< a M-k)-t{o))0(-k) <|0(-jfe)|, VkeR d , 

where ip G L 1 (R d , z). Thus, an application of the Lebesgue dominated con- 
vergence theorem yields 

lim/" dxe tA p(x)z(x) = [ z(dk)0(-k) lim e ^) d/2 (a(-fc)-a(o)) ) 



where 



lim e *(2^) d/2 (a(-fc)-a(o)) 



1, k = 

t"oc" 1 0, 7^ 

■ 

The next lemma shows that the existence of the arithmetic mean is stable 
under /^-convergence. In particular, it yields that the condition assumed in 
Proposition l5.5l on the Fourier transform of the intensity is sufficient to ensure 
the existence of the arithmetic mean. 

Lemma 5.11 Let z > be a bounded measurable function. Assume that 
there exist a total subset of L 1 {R d ,dx) and a C > such that for all p in 
that total subset we have 

lim R~ d / dxtp(x/R)z(x)=C / <p(x)dx. (31) 
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Then z has an arithmetic mean. In addition, equality l[37\) holds for C = 
mean(2;) and for any ip G L 1 (R d , dx). 

Proof. Let ip be an arbitrary function in L 1 (]R d ,<ix). Then for any e > 
there exists a if) in the aforementioned total set such that \\ifi — V'lU 1 < £ an d 
if) fulfills f l3T|) . By the following estimate, fl3Tl) also holds for <p: 



R d dxtp(x/R)z(x) — C / dxip(x) 



< 



dx (<f(x) — if>(x)) z(xR) 



\C\ 



dx(ip(x) — ip(x)) 



R 



-d 



dxip(x/ R)z(x) — C / dxif>(x 



< (|ML + |C|) y 



R- 



dxijj(x/R)z(x) — C / dxip(x 



Hence equality f[3"Tj) holds for any <p G L 1 (IR d ,(ix). Applying the result ob- 
tained till now to 



<Pr(x):=t r 1( ^ {r)) ' X f ^ 
0, otherwise 



r > 



(32) 



yields that the l.h.s. of (13TT) coincides with the mean(z). 



Corollary 5.12 Given a bounded measurable function z > the following 
two results hold: 

1. If z has arithmetic mean, then the limit in (3l\) exists for all ip G 
L 1 (R d , dx) and C = mean(z). 

2. If the Fourier transform of z is a signed measure, then z has arithmetic 
mean and 

mean(z) = (2vr)- d/2 z({0}). 



Proof. The first part is a direct consequence of Lemma fa. Ill using the total 
set of all ip r , defined as in (I32p . and their translates. 

If z is a signed measure, then due to Parseval's formula for all <p G 5(IR o! ) 
it holds 

R~ d i dx<p{x/R)z{x) = z{dk)<p(-kR) -»■ z({0})<p(0), R -> oo. 
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The limit exists because (p is continuous and decays quicker than any inverse 
polynomial. Hence, when R goes to oo, 0(—kR) converges to the charac- 
teristic function of the set {0}. The family (tp(—kR))R>o is dominated by 
an integrable function. The second part then follows by Lebesgue's domi- 
nated convergence theorem and an application of Lemma 15.111 for the total 
set S(R d ). M 



Remark 5.13 Let us underline the difference between £({0}) and the eval- 
uation of z as a function at zero. We explain this for the case when z is a 
bounded L 1 -function. The Fourier transform of z is then a continuous func- 
tion which we denote for the moment by z. Evaluation at zero makes sense 
in this case. However, we want to consider the Fourier transform as a gen- 
eralized function, 1.6. ? clS db linear form on a function space. We assume that 
this linear form is regular enough to be expressed by a signed measure. If, 
as in our case, z is an L 1 -function, then its Fourier transform is the measure 
z(dk) = z(k)dk. Thus £({0}) = 0, which does not necessarily coincides with 
*(°) = I** dxz{x). 

Below we list the Fourier transforms (interpreted as generalized func- 
tions, signed measures, respectively) of the examples given in Examples 15.71 
provided they exist and with the same enumeration: 

(i) (2vr) d / 2 z5 (dk). 

(iii) a continuous function for p = 1 or an L p ^ p_1 ' ) -function for 1 < p < 2 
multiplied in both cases by the Lebesgue measure. 

(iv) (27r) d / 2 (6o(dk) + ie/25 1 {dk) - ie/28-^dk)). 

(v) In the one dimensional case the Fourier transform is the following gen- 
eralized function z(k) = \^2n(z + z 1 )/25 (k) + i(z — z\) j 1 \ph\V{\jk\ 
where V(l/k) denotes the Cauchy principal value of 1/k. Using this ex- 
plicit formula the conclusion of Proposition 15.51 can be shown although 
the assumptions of Proposition 15.51 are not fulfilled. 

Applying the same technique as in Lemma 15.101 we can prove that the 
time asymptotic depends only on the behavior of z at infinity. 

Corollary 5.14 Let Z\,Zi > be two bounded measurable functions. If Z\ — 
z 2 G L 1 (IR d , dx), then the free Kawasaki dynamics with initial distribution n Zl 
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and the free Kawasaki dynamics with initial distribution tc Z2 have the same 
large time asymptotic limit. 

Proof. For all if G ViW 1 ) it follows from Lemma 15.101 and Corollary 15.121 
that 



dxe tA (e v - l)(x) Zl (x) - I dxe tA (e v -l)(x)z 2 (x) 
converges when t goes to oo to 



As we discussed in Example 15 .7\ the assumption on Z\ — z 2 implies that 
mean (24 — z 2 ) = 0. Hence 

lim / F* t (dry) e<*»> = lira / P* t (dj) e<*<* . 

By |Kal76t Theorem 4.2]), this is enough to show the required result. ■ 

Throughout this subsection, the study of the time asymptotic behavior 
of the free Kawasaki process with an initial distribution n z was based on the 
analysis of the Laplace transform 

^7r z (d 7 )E 7 [e^ x *>], (33) 

cf. Lemma 15.91 Actually, we had studied the time asymptotic behavior of 
the so-called empirical field corresponding to a <p G V(M. d ), 

n t (y?,X) := (<p, X$) = ^ 

xGX t 

The special role of the empirical fields is essential in the next subsection. 



5.2 Hydrodynamic limits 

In the sequel let z > be a bounded measurable function. In order to 
obtain a macroscopic description of our system, we rescale simultaneously 
the empirical field n t ((p, X) = ((p, X t ), tp G T>(M. d ), in space and in time. The 



32 



scale transformation in space is given by (ip, 7) — > e d (p(e-), 7), and in time by 
t — > e~ K t for some k > 0. To obtain non-trivial macroscopic density profiles, 
one has to scale the initial intensity as well, z — > z(e-). This scaling yields 
a scaling of the Laplace transform of the empirical field, in other words, the 
Laplace transform of the one-dimensional distribution of the scaled process. 
For each t > and each <p G T>(M. d ) we obtain from (1131) the following form 



*,(«.) (*y)e* fe E " KM (e"W 



1 ,7 



(34) 



exp fj^dx (e^ tA (e^-l)) 



(x)z(ex) ) . 



In the sequel we denote the scaled empirical field by 

4\cp,X) :=n e -« t (e d p>(£-),X)=s d (p(e-),X £ - H ) . 



(35) 



According to the independent movement of the particles, we are again able 
to reduce the study of the infinite particle system to an effective one par- 
ticle system. Again technical difficulties arise from the fact that the scale 
of the system size is much larger than the scale of space and time con- 
sidered in the empirical field. Actually, the system size is infinite. Under 
the additional assumption that the Fourier transform of the activity z is a 
signed measure, the hydrodynamic limit can be derived rather directly us- 
ing Fourier techniques, cf. Propositions 15.151 and 15.171 The general case of 
just bounded activities requires more technical involved considerations (post- 
poned to Proposition 15. 19j ) . 

Proposition 5.15 Let z > be a bounded measurable function such that its 
Fourier transform is a signed measure. For each t > the following limit 
exists for all p G ©(M^) 



lim / 7r*( e .)(d7) E 7 
e^o+ ' 



e ^ e) (^,x)" 


= / 




Jv 



5 Pt (du)e M 



(36) 



whenever one of the following conditions is fulfilled: 

afp := / dxXia(x) < 00, Vi = 1, 



,d, 
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and a/ 1 -* := (a[ , ■ • ■ , a^) ^ 0, then for k = 1 the limiting density pt is 
given (as a generalized function) by 



dx pt(x)(p(x) := / dxz{x + ta^ ] )ip(x), •p G T>(M. d ) 

2. 7/ a W = 0, and 

a ij ■— / dxXiXja(x) < oo, Vi, j = 1, • • • , d, 



then for k = 2 the limiting density p t is given (as a generalized function) 
by 



dx p t (x)ip(x) :- 



1 



(271)^/2 



dxz(x) / dke l{k ' x) e--z {ai2)k ' k) fi(k), ip G V(R d ) 



where denotes the d x d matrix with coefficients 



.(2) 



Proof. By ( 1341 . for each t > and each cp G X>(R d ), one has 
7r«( e .)(d7) E 7 



'r 

exp 



3 n< e) ( V ,X)' 



da;e £ " K ' A (e e ^ {£ - ) - l)(x) z(ea;) 



exp ^ tte e £_K ' A (e^^ - 1 - e VH) (*) *(e*) J 



•exp I / dxe d z(ex) ( e £ "^(^j ) (x 



(37) 



(38) 



Concerning ((37]), we observe that due to the contractivity property of the 



x z ex 



semigroup {e*~) t >a in L x (]R d ,(ix) (cf. Proposition I2.3P we have 
/ dxe £ - KtA (e ed ^ -l-£ d <p(e-) 
< \\e sd ^-l-eMe-)\\^ ld ,)\\z\\n 



(39) 



for more details see Lemma IA.1I in the Appendix |X] below. Thus, the proof 
reduces to check the existence of the limit of the exponential (|38p when e 
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converges to 0. For this purpose one shall apply the Parseval formula to 
the exponent in 0381) . The use of the explicit formula for the semigroup 

(e M ) t > 



' ! 1 - then yields 



z(dk) e £ 



^(27r) d / 2 (a(-efe)-a(0)) 



0(-k). 



(40) 



Here, observe that since for all k G M. d , Re(a(k) — a(0)) < (cf. Remark [23]) . 
one has 

;-H(2*r/*(a(-ak)-a(a))0(_k) < \0(~k) | , V k G R d , E > 0, 

where, in particular, (p G L 1 (M a! ,z). Therefore, due to the Lebesgue domi- 
nated convergence theorem, one may infer the existence of the required limit 
from the existence of the limit 



lim 



d(—ek) — a(0) 



(41) 



for each k G W 1 . 

Case 1: Let k — 1. Then, under the assumptions stated in 1, the function a 



is differentiate at 0, and thus (l4T|) exists with 

d(—ek) — a(0)) 



lim 



e ~ (2n) d / 2 

for each k G M. d . Hence the limit ( 138]) exists, being equal to 

11 1 ' J{k,(-x+taW)) r( _ 



z(dk) e it{k ' aW) 0(-k) 



dx z{x)(f{x — ta^ 1 ') 



0(-k) 



/ dxz(x + ta {1) )ip(x), peV( 



Therefore, the Laplace transform of the rescaled empirical field converges to 



exp I / dx z[x + ta^)(p(x) 1 
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This means that the limiting distribution is the Dirac measure with mass at 
z(x + 

Case 2: Now let k — 2. Then (f¥Tl) can be written as 



a(-ek)-a(0) 1 f , . e ^' fe > - 1 

ax a x 



e 2 (27r) rf / 2 ./ rod ~~ ~ v ~ y : 2 



1 



(27T) d / 2 



dot; CXj ( 3/ 



dec ci(cc^) I ^} ^ ^ ^ cCiCCjkikj ~\~ 

\ i,j=l 
( if o( £ 2 )\ 

2^ XiXjkikj + J 



where we have used the Taylor expansion of the function e*^'^ at x = and 
the assumptions for the Case 2. Since a G L^R , cfx), the latter converges 
to -l/2(2n)- d / 2 (a^k, k), and thus 



n< £ Vx)' 



lim / vr 2(£ .)(rf7)E 

e-*0+J T 

expf / 5(rfA;)e-5( a(2)fc . fc )^(-A;) 



exp 



-4tt7 / <Mrr) /" difee- i < fc ' a! >e-3< a(a)fc ' fc >^(-ife)^ . 

(27r) a! / 2 j R d V y 



Remark 5.16 For continuous differentiate z the limiting density pt{x) = 
z(x + ta^) obtained in Proposition 15.151 is the strong solution of the linear 
partial differential equation j^p t (x) = (a^\ Vp t (x)) = div(a^p t (x)) with 
the initial condition p = z - m the same way, the second case stated in 
Proposition 15.151 yields a limiting density which is the strong solution of the 
heat equation 

d , s if (2) 5 2 , x 
ij=i ■> 

with the same initial condition. 

Given the function a, we may decompose a into a sum of an even function 
p and an odd function q, a = p + q. Note that one always has p^ 1 ' = 0. 
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Submitting also the function a to a proper scale transformation (beyond the 
previous scale transformations in space and in time), one may complement 
the statement of Proposition 15.151 The limit (j4ip required in Proposition 
15. 151 suggests the scaling a e := p + eq and k — 2. Note that the assumptions 
on the function a (i.e., < a G L 1 (M d , dx)), carry over to a £ , e > 0. Hence all 
considerations done in Section [2] and following sections for the one particle 
operator A and its underlying dynamics still hold for the operator A e defined 
by (jl]) with a replaced by a E . Denote by 

:= / dxxa(x)=e~ 1 / dxxa £ (x) = / dxxq(x). 



Proposition 5.17 ("weak asymmetry") Let z > be a bounded measurable 
function such that its Fourier transform is a signed measure. Under the above 
conditions, if ^ = G M d and pf) < oo for every i,j = l,---,d, 
then for each t > and each G V(R d ) the following limit exists, and it is 
given by 



lim / 7r, (e .)(d7)E 7 [ e ^ } ^ x )l 

1 f dxz(x) I dke^e-^^-^^^ik) 



cxp 



(2vr) d / 2 



Proof. Since all statements for the operator A and the corresponding 
process also hold for A £ , in particular, one finds that (e tAe ) t >o is also a 
contraction semigroup on L 1 (M d ,dx). Similar arguments as in the proof of 
Proposition 15.151 reduce the proof to the analysis of existence of the limit 

^ S(-rt) -g(0) = fl-E*) -flO) + , im fl-E*)- j(0) fceR , 



Remark 5.18 Similarly to Remark 15.161 one may then conclude that for 
continuous differentiable z Proposition 15.171 leads to a limiting density p t 
which is a strong solution of the partial differential equation 



d , , , m , ^ lv^ (2) 9 2 



1 d 

-p t[ x) = div(a^ Pt (x)) + g E ^'d^P^ 

i,j=i ■> 



with the initial condition P q = z. 
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Proposition 5.19 Assume that a has all moments finite. Then the results 
stated in Propositions \5.15\ and \5.17\ hold for all non-negative bounded mea- 
surable functions z. 



Proof. We note that in the first part of the proof of Propositions 15.151 and 
15.171 we have only used the boundedness of z. Therefore, it remains to prove 
the convergence of the exponent in (138]) . namely, 



dxz(x) (e £ KtAe <f(e 

jR d ^ 

Let if) G L 1 (R d , dx) be given. Then 

dxz(x) (e £ ~ KtAe (p(e 
< e a I dx \z(ex 



(42) 



< \\4uS C 



^£~ K tA e 



dxz(x) (e £ RtAe if)(e-\ 

/ Vfc/ jR d ^ 

^(eO-VCe-)) 



L 1 (R d ,dx) 



As (e tA z )t>o is a L 1 (R d , (ix)-contraction semigroup, the latter can be bounded 
by 



< z\Le c 



dx \<f(ex) - ip{ex)\ = \\z\\ u \\(p - ^IU^rc^)- 



Therefore, it is enough to consider (j42p for (p from a total subset of L 1 (M d , dx). 
Let us consider the set of all functions in L 1 (R d , dx) with Fourier transform in 
V(R d ). This set is total in L 1 (R d , dx), because V(R d ) is dense in S(R d ), the 
Fourier transform is continuous in iS(M d ), and iS(lR d ) is dense in L 1 (M' i , dx). 
Let ip be such a function. In order to prove the convergence of (l4"2l) it is 
enough to show the convergence of (e £ KtAe (f(£-))(j) in L 1 (R d ,dx). For this 
purpose it is sufficient to show the following convergence of the Fourier trans- 
form of the latter expression, 



-t(27r) d / 2 (d £ (~ek)~a(0))<n(- 



<p(-k) 



e ft{.« l *)-^(.P)M)^_ Jfe)i £ ^ (43) 



in <S(R d ). The exponent in the l.h.s. of f|4*3|) may be written as 



a £ (-ek) - a(0) = -e(Va e (0), k) 



(27t) d / 2 



'1-s) 



(k,x) 2 e is£ ^ x) a £ (x)dxds.(U) 
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The first summand on the right hand side is of order e K , because e(Va e (0), k) = 
j^Xdjz J^d dx(k, x)q(x). Hence all derivatives of this expression grow at most 
polynomially and all derivatives of third order and bigger are of order e or 
smaller. Hence only derivatives of the first and second order have to be 
computed carefully in order to see that 

e ~£ 2 - K Jo(l-») I ud {k,x)*e™e< k >*)a s {x)dxds _ e -<?Zl± {a (*) k)k) ^ 

are polynomially bounded of order e. Summarizing, (143p and all its deriva- 
tives converge locally uniformly. As (p G V(M. d ), ( l4~3i) also converges in iS(lR d ). 
Appendix [B] provides the tools for a more detailed calculation. ■ 



6 Non-equilibrium dynamics 

In this section we widen the class of initial distributions to measures far 
from equilibrium, that is we consider all probability measures /i on as ini- 
tial distributions subject only to a mild mixing condition. This means that 
we consider the processes constructed as in Section [3] but not necessarily 
with a Poissonian initial distribution. Assuming enough mixing of the initial 
measure fi, namely (1471) . we are able to generalize, incorporating ideas from 
[DSS82J, Proposition 15.51 in Subsection 16.11 and Proposition 15.191 in Subsec- 
tion o 

We formulate the mixinig requirement in terms of the second Ursel func- 
tion (factorial cumulant), which can be expressed in terms of the first and 
second correlation function (factorial moments) defined in Subsection 13. 2\ 
namely 

uf\x,y) := k„({x,y}) - k^({x})k^({y}). 

We denote in the following the first correlation function x \— > k^{{x}) 
by pp. The condition on the second Ursel function, (|47p . is a rather weak 
mixing or decay of correlation condition. In Subsection 16.31 we show that this 
condition is fulfilled, in particular, by Gibbs measures in the high temperature 
regime. It will also hold beyond that regime, cf. e.g. [BMPK04J. 



6.1 Large time asymptotic 

Recalling (j!7p . the Laplace transform of the one- dimensional distribution P^ t 
can be expressed in terms of an one particle system, i.e., for all non- negative 
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/ e S(R d ) 



e B {e tA {e- f -l), 1 )ii{d 1 ) 



e 



3 {ln(e tA (e--f~l)+l), 7 ) 







//(d7). 
(46) 



Proposition 6.1 Lei n be a measure on which has first and second cor- 
relation function. Assume that \x fulfills the following mixing condition 



sup x&K d I ujjp(x,y)dy<oa. 



(47) 



In addition, we assume that the Fourier transform of the first correlation 
function p M is a signed measure. Then, the one- dimensional distribution P* t 
converges weakly to 7r mean ( p \ when t tends to infinity. 

Proof. Given an non-negative / G <S(IR d ) such that — 1 < e~f — 1 < 0, let 
ip := 1 — e~*. Using ( 146]) and \e~ x — e~ y \ < \x — y\ for x, y > one obtains 
that 



,<ln(e tA (e--f-l)+l),7)^ c ^ _ e moan(p M ) J Rd (e~i ' -l)(x)dx 







< 



(ln(— e (f + 1),7) - mean(p At ) / tp(x)dx 



As < x — ln(x + 1) < x 2 for — 1/2 < x < 0, ||e* A (/9|| u tends to zero for 
t — > oo, because G L 1 (IR c( , dx) and 



e*\>(*)l < 



(2vr) rf / 2 



d j fce W , ^B*(a(*)-a(o))|0( A .)| > 



and (e* A ) t >o is an L 1 (M d , (ix)-contraction we get that 

f |(ln(e^ + l), 7 )-(e i V7)|Mrf7)< f (e tA <p(x)) 2 P^)dx 
Je JR d 

< \\e tA ip\\ u e tA (p(x)p fi (x)dx < ||e ty V|| u ||y?|| L i ||p M || u . 



This means it is sufficient to estimate 



/ (e tA ip, 7) - mean(p M ) / <p{x)dx 
Je Jwi d 
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< 



(e 4 V7>- / e tA cp(x)pJx)dx 



tA 



+ 



e tA (p(x)p fJi (x)dx — mean(p M ) / (p(x)dx 



The last term converge to zero because of Lemma 15.101 and Corollary 15.121 
Due to the decay properties of the covariance function ff4T|) 



e 



(e>,7)- e tA if(x)p ll (x)dx 



< 



i tA (p(x)e tA (p(y)u^{x,y)dxdy + / {e tA ^{x)f Pfl (x)dx 



< He^Vye^llx! ( sup / uW(x,y)dy+\\ Plt \\ 



tA, 



,( 2 )/ 



which implies the result, as e is an L 1 (IR , dx) contraction and ||e V^IU 
converges to zero. ■ 



6.2 Hydrodynamic limits 

As in Subsection 15. 21 we want to study the rescaled empirical field n[ e ^ (<p, X) = 
e d (ip(e-), X £ -«t). Since we do not have any longer a natural parameter as- 
sociated to the initial measure, one cannot formulate something like slowly 
varying intensities. However, one sees that a possible framework is to work 
with a quite arbitrary sequence of initial measures (p £ ) £> Q. The main restric- 
tion on this sequence is that one has to assume a particular convergence for 
the first correlation measure described below in more details and the mixing 
condtion (I47p uniformly in e. In Corollary 16.41 we prove that these conditions 
are fulfilled by Gibbs measures in the high temperature regime for slowly 
varying intensity z(e-). Furthermore, the limit is identified. 

Theorem 6.2 Assume that a has all moments finite. Let (p e ) £> o be a se- 
quence of measures on 6 such that p /Ue is uniformly bounded in e, the limit 
lim £ ^ + Pn e ({x/e}) =: P o(x) exists for all x G M d and the following mixing 
condition 

sup / u ( £(x,y)dy < oo 

x£R d ,e>0 JR d 
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holds. Then, for each t > 0, the following limit exists for all non-negative 

cp e v(R d ) 



e->0-t 



lim / p s (dj)E^ 



S pt (du)e 



-(tp,uj) 



(48) 



v 



whenever one of the following conditions is fulfilled: 
,(D _ ... „(U 



a\ a d ) 7^ 0, then for k = 1 t/ie /zmzt ffllfy holds for 



1. If a 

p t equal to the strong solution of J^p t (x) = div(a^ p t (x)) with initial 
condition po- 

2. If = 0, then for k = 2 the limit HM holds for p t equal to the strong 
solution of ^p t (x) = | Yutj=i a ij eSdx Pt( x ) initial condition p . 



3. As in Proposition \5.11\ let p and q be the even and the odd part of a. 

Then, for 



(!) = 



Define a e = p + eq and assume that ^ q 
k — 2, the limit $B\) for the dynamics w.r.t. a e holds for p t equal to 
the strong solution of §- t p t (x) = div(a (1) pt(») + \ J2i,j=i a if aS^ pt ^ 
with initial condition p . 

Proof. Given an non-negative <p G 5(IR d ) according to ( l4"6j) one may write 
flUD as 



p £ (d-f) e 



<ln(e^- K ^(e- d ^-)_l) + l), 7 ) _ kd v{x)pt{x)dx 



e 



(49) 



Using |e x — e y \ < \x — y\ for x, y > one can bound this by 



(In j>-"t^( e -*W0 _ l) + l) j7 ) + / v (x)p t {x)d. 



e 



Let us show that it is sufficient to consider 



fi e {dj) 



{e £ K ' Ae eVH,7>- / <p{x)p t {x)dx 
indeed, proceeding as in the proof of Proposition 16.11 we get 

J | (In (>-^( e -<VH _ i) + i^ )7 ) + ( e ^eMe-)n)\^(dl) 



(50) 
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< 



f In ( e £ ~* tA * (e- £d ^) - 1) (x) + l) - e £ ~ K 



l)(x) p^{ dx ) 



+ 



l){x) + e*- KtA *e d ip{e-){x) p, £ (dx) 



< 



e £ ~ KtA > I e 



-^('0-l+eVeO)|(aOftJ<&0 



< 



1) 



3 e K tA e / -£ d ¥5(£-) 



L 1 



+ 



l+e d V {e-)) J/vJ u . 



L 1 



The latter expression is at least of order e d according to Lemma IA.1I 
Thus it remains to consider (I50|) which can be bounded by 



e 



( e £ ^e<VH,7> - I (e £ ~ KtA ^(e-) ) {x/e)p, e {x/e)dx 



,e- K tA E 



ip(£-) ) (x/£)p fls (x/e)dx - / ip(x)p t (x)dx 



Pe(dj) 

(51) 



In Proposition 15. 19l we show that e £ * e ip(e-)(x/s) converges in iS(IR d ) to 

By assumption, P/j, s ({ x / £ }) converges, in particular, in || • || 0,-^-1,2 5 cf- (|62|) - 
and thus one obtains 



lim / e 



; HA ^{e-){x/e)p lie {{x/e})dx 



Hence the second summand in (1511) converges to zero. The first summand can 
be bounded by the second Ursel functions as in the proof of Proposition 16.11 



e 



< 



e £ ~ KtA *e d ip(e-) 



,e- K tA E 



<p(e-)) (x/e)p^(x/e)dx 



Pe(dj) 



e £ - KtA °e d <p{e-] 



L 1 



sup 



x£R d .E>0 



u ( ^(x,y)dy+ \\p^\\ u 
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6.3 Application to Gibbs measures 

In this subsection we prove that the hypothesis for the results of the previ- 
ous subjection are fulfilled for a concrete class of non-equilibrium measures, 
namely, for Gibbs measures in the high temperature low activity regime. 

In order to recall the definition of a Gibbs measure, first we have to 
introduce a pair potential V : R d -> R U {+00}, that is, a measurable 
function such that V(-x) = V{x) E R for all x E R d \ {0}. For 7 E T 
and x E R d \ 7 we define a relative energy of interaction between a particle 
located at x and the configuration 7 by 

if ^2\V(x-y)\ < 00 

ye-y 

otherwise 

A probability measure /i on T is called a Gibbs measure corresponding to 
V, an intensity function z > 0, and an inverse of temperature j3 whenever it 
fulfills the Georgii-Nguyen-Zessin equation |NZ79| Theorem 2] 

I li{d-f) Y^H(x,i) = f nidrf) I dxz(x)H(x,-fU{x})e- mx ^ (52) 
Jr x£7 Jr Jr* 

for all positive measurable functions H : R d x T — ► R. This definition is equiv- 
alent to the definition via DLR-equation, see [Geo76t INZ791 IPuz81t IKun99j . 
We observe that for V = (1521 reduces to the Mecke identity, which yields 
an equivalent definition of the Poisson measure ir z |Mec67l Theorem 3.1]. We 
also note that for either V = and z not being a constant or V 7^ 0, a Gibbs 
measure neither is a reversible nor an invariant initial distribution for the 
free Kawasaki dynamics under consideration. In order to have thermody- 
namical behavior we assume that V is stable, i.e., there exists a B > such 
that ^2^ x y y Cv V(x — y) > — B\q\ for all configurations r\ E T . Furthermore, 
we shall assume that the parameters /3, z are small (high temperature low 
activity regime), i.e., 

\\z\\ u e^ B+1 C((3) < 1, 

where C((3) := J Rd dx\e~ l3V ^ — 1|. These conditions are, in particular, suf- 
ficient to insure the existence of Gibbs measures, cf. |Rue69] . Moreover, the 
correlation functions corresponding to such measures exist and fulfill a Ruelle 



E 

2/67 



V(x-y), 



+00, 
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bound denned in Section [3J and thus, as noted there, they are supported on 
6. 

In the high temperature low activity regime one has rather detailed in- 
formation about the Ursell functions (factorial cumulants) : Tq —>■ E 
corresponding to /i, which are bounded measurable functions (also called 
Ursell functions), i.e., for all / G X>(lR d ) holds 

I e<^V(d 7 ) = exp ( / T[(e f W - l)u IM (r])X(dr]) \ . (53) 
Jr V r ° yen J 

The function x —>■ m m ({x}) coincides with the first correlation function of \x. 

We present the results necessary for the following. For further details, see 
e.g. jDSI75j . |MM91j and see also [KunOlj . 

The Ursell functions can be expressed in terms of a sum over all connected 
graphs weighted by the Meyer-functions 

fc ® : = E II {e-M*-y)-i) (54) 
where G c (0 denotes the set of all connected graphs with vertex set £: 

Actually, the following bound is the key result of the cluster expansion of 
Penrose-Ruelle type 

\ m <e^\ II {e~ mx - y) -ll (56) 

T&T{0 {x,y}£G 

where T(£) denotes the set of all trees with set of vertices £. This leads to 
the following integrability bound 

/ \u^({x,yi, • • • ,y n })\z(y 1 )dy 1 . . . z(y n )dy n 

JR dn 

< e ^ B+1 ^ (\\z\\ u C{(5)f Y (n + m , + 1)! (e^ B+l \\z\\ u C((3)) m (57) 

/ — ' ml 

m=0 

In particular the mixing condition (T4T1) of Proposition 16.11 and Theorem 16.21 
holds. 

We show that Gibbs measures in the high temperature regime with a 
translation invariant potential fulfill the assumptions of Proposition 16.11 
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Corollary 6.3 Let z > be a bounded measurable function which Fourier 
transform is a bounded signed measure. Let p be a Gibbs measure correspond- 
ing to a translation invariant potential V described above, inverse temperature 
j3 and activity z which are in the high temperature low activity regime. Then 
the first correlation function p M has as Fourier transform a measure and the 
arithmetic mean 

1 °° 1 f ~ 

mean(p M ) = , V -r / k r {pi, . . .,p n )z{{pi+. . .+p n })z(dp 1 )-. . .-z(dp n ), 

where 

k r (y u ...,y n ) :=£ J] ( e -W*-*») - 1). 
Geg c {xi,x2}eG 

i/ere (? c denotes the set of all connected graphs with vertex set (0, y±, . . . , y n ). 
As a consequence, all assumptions of Proposition WA\ are fulfilled. 

Proof. Due to the translation invariance of V and cluster expansion one 
can rewrite the first correlation function as 

p li (x) = ^2— Yl II {e~ pv{vi ~ y2) -l)z{x l -x)-...-z{x n -x)z(x)dx v ...-dx r 

where Q c denotes the set of all connected graphs with vertex set (0, yi, . . . , y n ) := 
(0, Xi — x, . . . , x n — x). By fl56l) and fl57|) the function k r is integrable, and 
thus k r is a continuous function which decays to zero at infinity. Hence 

k r {x x -x,...,x n -x) = (2vr)- nd / 2 f dp x ... dp n dpe ipixi e ipnXn e ipx 

JK dn 

■k r (p 1 ,...,p n )S(p-p 1 - ... -p n ). 

To compute the Fourier transform of p^ in the weak sense it remains to 
compute the Fourier transform of 

ip(x)z(xi — x) ■ . . . ■ z{x n — x)z(x) 
= (2rr)- (n+2)d/2 [ z(d Pl )e tpiXl z{dp n )e iPnXn dp e ipx 

JM. dn 

e -i( Pl+ ... +Pn)x ^ 



46 



for any ip G D(M. d ). Summarizing, we obtain that 

/ dxip(x)p fl (x) = ^ —i I K d Pi) ■ ■■■■ z(dp n )z(dp) 

JR d n=Q n - JR dn 



■fc r (pi, • • • ,Pn)0{Pl + • • • +Pn -p)- 



As in Subsection 15.21 we consider initial measures with a slowly varying 
intensity, i.e., Gibbs measures corresponding to /3, V, and z(s-), which we 
denote by \x e . As \z\ u is unchanged, all scaled measures \i e remain in the high 
temperature low activity regime and the bound (15Tj) holds uniformly in e > 0. 
For Gibbs measures which are not Poisson measures, the first correlation 
function is not any longer just the intensity. The function appearing as initial 
value in the limiting partial differential equation is the scaling limit of the 
first correlation function and not just the unsealed activity. Let us describe 
what is the scaling limit of the first correlation function. Denote by p^ qm 
the correlation function corresponding to the Gibbs measure with constant 
activity c, temperature (3 and potential V. Due to the translation invariance 
of V this correlation function is a constant function. Given a function z > 
denote by x i— > the function which associates to each x the constant 
value of the first correlation function of the Gibbs measure with constant 
activity z(x). This function is the scaling limit of the first correlation of \i e . 
Note that p M (x) = u^({x}). 

Corollary 6.4 Given a bounded measurable function z > 0, a potential V , 
and an inverse temperature (3 fulfilling the conditions of the high temperature 
and low activity regime, the corresponding Gibbs measures fulfill all assump- 
tions of Theorem \6.2 . Moreover, 

Pa(x):= lim . Ufl£ ({x/e}) = p e Tl (58) 

£— >0+ v ' 

Proof. According to the cluster expansion of and the translation invari- 
ance of V we obtain that 

oo „ n 

n=l ' jRd " Gee c ({0,...,n}) {i,j}eG 1=1 

(59) 

where x/q := 0. Due to fl57|) the above expression is uniformly integrable in e 
and in x, yielding the claimed result. ■ 
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A Estimates for hydrodynamic limits 

Lemma A.l Let if G S(R d ) and < e < 1 be given. Then 

|| e u (e ^H _ 1)11^^^ < IMUi^e 11 ^ 1 ^) (60) 

\ e e^e-) _ 1 - £ V(£-))IUHM^) < ^1^lliu^<V V " il(Md ' da;) 



Proof. On the one hand (e tA ) t >o is a contraction semigroup with respect 
to the the supremum norm. So 

00 nd 

\\e tA (e^ - l)|| u < \\e^ - l\\ u < V ^HMI" 

z — ' n! 

n=l 

On the other hand (e tA ) t >o is a contraction semigroup on L 1 (M. d , dx). There- 
fore, the left-hand side of (!60l) is bounded by 

f I^V*) _ 1| e d(n-l) 

/ dx Td IMUWZ^— Za~ WvWlHdx) < ll^lliH^) e 

The last inequality follows by a similar computation. ■ 



B Norm estimates in S 



)d 



Let us introduce the following two equivalent systems of norms for the locally 
convex topological vector space iS(IR d ) for AgM and M > 0: 



A,M,u 



AM3 



swp\D«f(x)\(l + \x\ 2 ) M (61) 

\a\<A 

£ (l d \Dy(x)\ 2 (l + \x\>) M d x y 2 (62) 



\a\<A 



Lemma B.l Let f u f 2 G C°°(R <i ; C) 6e two C°° -functions with non-positive 
real part such that for an A G N and a M > one /ias ||/i|U,-M,u < oo. 
T/ien t/iere exists a constant C depending on A and on \\f%\\A-M,u such that 

\\e h - e /2 || A ,-(A+i)M,2 < C\\fi - / 2 |U,-m,2- (63) 
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Proof. By induction one finds 

£>V-1)(*)=E E f[D ai f(x)(e f to-l). (64) 

k=1 ai,...,a fe gN d J =1 
a\+...+a k =ct 

Thus, using telescope sums one can write the following difference as 

\a\ k (l-l k 

D»(e* - e*)(x) = E E E [\{D^ h{x)D^ h{x) J] ^/ 2 (x) 

ai+...+o fe =o 
l-l k \ 

j=l j=Z+l / 

k 

+ IJ D ai h{x){e Mx) - e f2(x) ). 

3=1 

Using Cauchy-Schwarz inequality one may estimate this by 

\ 1/2 



\D«(e^ - e^)(x)\ 2 (1 + \x\Y (A+1)M dxj 



\a\ k 



< EE \\m\aUu\\m\X-mM - h\\A,-M,2\\e h \ 



k=l 1=1 



+ll/2|ll-M,J|e /l -e /2 |U,_M, 2 

< iar(i+ii/iiu-M, u +ii/ 2 iu-M, u ) |a| 

(Wfi -/2|U,-M,2||e /l ||u+ \\e fl -e /2 || ,-M, 2 ) • 

The non-positivity of the real part of fi and the properties of the exponential 
function yield 

||e* -e A 110,-^,2 < ||/i-/2||o,-m,2. (65) 

Putting the additional sum out by triangle inequality and using non-negativity 
one obtains 

\\e h - e h \\ A ,-{A+i)M,2 (66) 

< A 3 (l + ||/ 1 |U,-Af,u+ Wf2\\A,-M,u) A ||/l " / 2 |U,-M,2. (67) 
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Lemma B.2 Let f, f\, f 2 G C°°(IR d ; C) be such that fx, f 2 have non-positive 
real part and for an A G N and a M > one has \\fi\\A-M,u < oo and 
for a M G R one /ias ||/|U,m ,« < °°- Then for any A G N and /or any 
< M' < Mq — (A + 1)M there exists a constant C depending monotonically 
on A, \\f\\ A ,Mo,u, and on \\fi\\A,-M,u such that 

|| /e /i _ f e J>\\ A)M , >2 < CWfr - f 2 \\ A ,~M,2. (68) 

Proof. Due to triangle inequality and product rule one obtains 

\D a (f(e fl -e /2 ))| < \D ai f D a *(e fl -e h )\. 



ai+ct2— ot 



Then using triangle inequality and the supremums norm one finds 

l2\M -(A+l)M 



< 



£ (| |D-/(x)D-(e^-e^)(x)| 2 (l + |x| 2 ) M °- (A+1)M ) 



1/2 



II .ft .ft II 

A,-(A+1)M,2- 



< A||/|| AjMoiU ||e^-e^ 



Acknowledgment 

We would like to thanks N. Jacob, E. Lytvynov, M. Rockner, W. Hoh for 
fruitful discussions, T. Hilberdink for mentioning the theory slowly varying 
functions and T. Funaki for making us aware of [DSS82j . Yu. K., T. K., 
M. J. O. thanks CCM (University of Madeira) for the pleasant ongoing hos- 
pitality. Financial support from DFG through the SFB 701 (Bielefeld Univer- 
sity) and FCT through, POCI-2010, FEDER, PTDC, PTDC/MAT/67965/2006 
are gratefully acknowledged. 



References 

[AKR98] S. Albeverio, Yu. G. Kondratiev, and M. Rockner. Analysis and 
geometry on configuration spaces. J. Fund. Anal, 154(2) :444- 
500, 1998. 



50 



[BCC02] L. Bertini, N. Cancrini, and F. Cesi. The spectral gap for 
a Glauber-type dynamics in a continuous gas. Ann. Inst. H. 
Poincare Probab. Statist, 38:91-108, 2002. 

[BK88] Yu. M. Berezansky and Yu. G. Kondratiev. Spectral Methods 
in Infinite- Dimensional Analysis. Naukova Dumka, Kiev, 1988. 
(in Russian). English translation, Kluwer Academic Publishers, 
Dordrecht, 1995. 

[BMPK04] F. Baffioni, I. Merola, E. Presutti, and T. Kuna. A relativized Do- 
brushin uniqueness condition and applications to Pirogov-Sinai 
models. Preprint 04-107, mp-arc, 2004. 

[Bog46] N. N. Bogoliubov. Problems of a dynamical theory in statistical 
physics. Gostekhisdat, Moskau, 1946. In Russian. English trans- 
lation in J. de Boer and G. E. Uhlenbeck, editors, Studies in 
Statistical Mechanics, volume 1, pages 1-118, Amsterdam, 1962. 
North-Holland. 

[DMP91] A. De Masi and E. Presutti. Mathematical methods for hydro- 
dynamic limits, volume 1501 of Lecture Notes in Mathematics. 
Springer- Verlag, Berlin, 1991. 

R. L. Dobrushin. On Poisson's law for distribution of particles 
in space. Ukrain. Mat. Zh., pages 127-134, 1956. (In Russian). 

J. L. Doob. Stochastic Processes. John Wiley & Sons, New York 
and London, 1953. 

M. Duneau, B. Souillard, and D. Iagolnitzer. Decay of correla- 
tions for infinite-range interactions. J. Math. Phys, 16 (8): 1662- 
1666, 1975. 

R. L. Dobrushin and Ra. Siegmund-Schultze. The hydrody- 
namic limit for systems of particles with independent evolution. 
Math. Nachr., 105:199-224, 1982. 

S. N. Ethier and T. G. Kurtz. Markov Processes. Characteri- 
zation and Convergence. Wiley Series in Probability and Math- 
ematical Statistics. John Wiley & Sons, NewYork, Chichester, 
Brisbane, Toronto, and Singapore, 1986. 



[Dob56] 
[Doo53] 
[DSI75] 

[DSS82] 

[EK86] 



51 



[Fri87] J. Fritz. Gradient dynamics of infinite point systems. Ann. 
Probab., 15(2):478-514, 1987. 

[Geo76] H. 0. Georgii. Canonical and grand canonical Gibbs states for 
continuum systems. Comm. Math. Phys., 48:31-51, 1976. 

[G1681] E. Glotzl. Time reversible and Gibbsian point processes. I. 

Markovian spatial birth and death processes on a general phase 
space. Math. Nachr., 102:217-222, 1981. 

[JacOl] N. Jacob. Pseudo-Differential Operators and Markov Processes, 
volume I. Fourier Analysis and Semigroups. Imperial College 
Press, London, 2001. 

[Kal76] O. Kallenberg. Random Measures. Academic Press, London, New 
York, second edition, 1976. 

[KK02] Yu. G. Kondratiev and T. Kuna. Harmonic analysis on config- 
uration space I. General theory. Infin. Dimens. Anal. Quantum 
Probab. Relat. Top., 5(2):201-233, 2002. 

[KKK04] Yu. G. Kondratiev, T. Kuna, and O. Kutoviy. On relations be- 
tween a priori bounds for measures on configuration spaces. In- 
fin. Dimens. Anal. Quantum Probab. Relat. Top., 7(2): 195-213, 
2004. 

[KKO02] Yu. G. Kondratiev, T. Kuna, and M. J. Oliveira. Analytic aspects 
of Poissonian white noise analysis. Methods Funct. Anal. Topol- 
ogy, 8(4):15-48, 2002. 

[KKO06] Yu. G. Kondratiev, T. Kuna, and M. J. Oliveira. Holomorphic 
Bogoliubov functionals for interacting particle systems in contin- 
uum. J. Funct. Anal, 238(2):375-404, 2006. 

[KKZ06] Yu. Kondratiev, O. Kutoviy, and E. Zhizhina. Non-equilibrium 
Glauber type dynamics in continuum. J. Math. Phys., 
47(11):113501, 2006. 

[KL05] Yu. Kondratiev and E. Lytvynov. Glauber dynamics of contin- 
uous particle systems. Ann. Inst. H. Poincare Probab. Statist., 
41(4):685-702, 2005. 



52 



[KLR03] Yu. Kondratiev, E. Lytvynov, and M. Rockner. The heat semi- 
group on configuration spaces. Publ. Res. Inst. Math. Sci., 
39(l):l-48, 2003. 

[KLR07] Yu. G. Kondratiev, E. Lytvynov, and M. Rockner. Equilibrium 
Kawasaki dynamics of continuous particle systems. Infinite Di- 
mensional Analysis, Quantum Probability and Related Topics, 
7(4): 185-209, 2007. 

[KLR09] Yu. G. Kondratiev, E. Lytvynov, and M. Rockner. Non- 
equilibrium stochastic dynamics in continuum: the free case. 
Condensed Matter Physics, 11(4(56)):701-721, 2009. 

[KSSU98] Yu. G. Kondratiev, J. L. Silva, L. Streit, and G. F. Us. Analysis 
on Poisson and Gamma spaces. Infin. Dimens. Anal. Quantum 
Probab. Relat. Top., 1(1):91-117, 1998. 

[Kun99] T. Kuna. Studies in Configuration Space Analysis and Appli- 
cations. PhD thesis, Bonner Mathematische Schriften Nr. 324, 
University of Bonn, 1999. 

[KunOl] T. Kuna. Properties of marked Gibbs measures at high tem- 
perature regime. Methods Fund. Anal, and Topology, 7:33-54, 
2001. 

[Kun09] T. Kuna. Bochner's theorem for point processes via Bogoliubov 
functionals. In preparation, 2009. 

[Len75a] A. Lenard. States of classical statistical mechanical systems of 
infinitely many particles I. Arch. Rational Mech. Anal, 59:219- 
239, 1975. 

[Len75b] A. Lenard. States of classical statistical mechanical systems of 
infinitely many particles II. Arch. Rational Mech. Anal, 59:241- 
256, 1975. 

[Mec67] J. Mecke. Stationare zufallige MaBe auf lokalkompakten 
Abelschen Gruppen. Z. Wahrsch. verw. Gebiete, 9:36-58, 1967. 

[MM91] V. A. Malyshev and R. A. Minlos. Gibbs Random Fields: Cluster 
expansions. Kluwer Academic Publishers, Dordrecht, Boston, 
and London, 1991. 



53 



[NZ76] X. X. Nguyen and H. Zessin. Martin-Dynkin boundary of mixed 
Poisson processes. Z. Wahrsch. verw. Gebiete, 37:191-200, 1976. 

[NZ79] X. X. Nguyen and H. Zessin. Integral and differential character- 
izations of the Gibbs process. Math. Nachr., 88:105-115, 1979. 

[Puz81] J. Puzicha. Existenz und Charakterisierung von Gibbsmafien. 
Master thesis, University of Bielefeld, 1981. 

[Rue69] D. Ruelle. Statistical Mechanics. Rigorous Results. Benjamin, 
New York and Amsterdam, 1969. 

[Rue 70] D. Ruelle. Superstable interactions in classical statistical me- 
chanics. Comm. Math. Phys., 18:127-159, 1970. 

[Sin89] Ya. G. Sinai, editor. Ergodic Theory with Applications to Di- 
namical Systems and Statistical Mechanics, volume II of Ency- 
clopaedia Math. Sci., chapter III. Dynamical systems of statis- 
tical Mechanics and Kinetic Equations. Springer- Verlag, Berlin 
Heidelberg, 1989. 

[Spo86] H. Spohn. Equilibrum fluctuations for interacting Brownian par- 
ticles. Comm. Math. Phys., 103:1-33, 1986. 



54 



